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ABSTRACT 


The  effects  of  curvature  and  the  resulting  centrifugal  instabilities  on  local  heat  transfer 
distributions  are  studied  in  a  curved  channel  at  Dean  numbers  ranging  from  50  to  200.  The 
channel  has  a  rectangular  cross  section  of  1.27  cm  by  50.1  cm  giving  an  aspect  ratio  of  40 
to  1.  Flow  is  heated  in  a  straight  portion  of  the  channel  prior  to  the  curved  portion  in  ordci 
to  obtain  flow  which  is  hydrodynamically  and  thermally  fully  developed.  All  baseline  tests 
confirm  techniques  employed  and  qualify  flow  behavior.  These  consist  of  energy  balance 
checks  and  comparison  of  results  from  the  straight  section  to  numerical  and  analytic 
solutions.  Nusselt  numbers  in  the  curved  section  initially  show  an  abrupt  decrease  after  the 
imposition  of  the  stabilizing  influences  of  convex  curvature.  These  are  followed  by  a 
gradual  increase  as  centrifugal  instabilities  and  Dean  vortices  form  and  develope.  Spacially 
resolved  results  also  show  significant  surface  Nusselt  number  variations  across  the  span  of 
a  vortex  pair,  especially  on  the  concave  surface.  On  the  convex  surface,  local  Nusselt 
numbers  are  much  more  apt  to  be  spanwise  uniform. 
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I.  INTRODUCTION 


A .  DESCRIPTION  OF  THE  PROBLEM 

The  effects  of  curvature  and  the  resulting  centrifugal  instabilities  on  local  heat 

transfer  distributions  are  studied  using  a  curved  channel  with  an  internal  aspect  ratio  of 
40  to  1.  Results  are  presented  for  Dean  numbers  ranging  from  50  to  200,  resulting  in 
laminar  flow  which  sometimes  contains  pairs  of  Dean  vortices.  Other  work  at  the  Naval 
Postgraduate  School  in  this  area  addresses  the  development  of  vortex  flows  due  to 
centrifugal  instabilities  and  the  effects  of  these  secondary  flows  on  transition  from 
laminar  to  turbulent  behavior  without  considering  heat  transfer.  Ligrani  and  Niver  [Ref. 
1]  describe  some  results  from  this  study  which  are  also  obtained  in  a  curved  channel 
with  a  similar  aspect  ratio  of  40  to  1 .  These  authors  also  review  other  non-therma! 
studies. 

Several  other  studies  focus  on  heat  transfer  in  channels  with  turbulent  flow.  Johnson 
and  Launder  [Ref.  2]  present  local  heat  transfer  results  from  a  channel  with  a  180 
degree  bend  and  square  cross  section.  They  emphazise  that  instabilities  from  curvature 
produce  very  complex  examples  of  turbulent  flows  due  to  substantial  secondary  flows. 
Additional  studies  are  described  by  Chang,  Humphrey,  Johnson  and  Launder  [Ref.  3],  who 
employ  a  channel  with  a  180  degree  bend  of  square  cross  section,  and  by  Brinch  and 
Graham  [Ref.  4],  who  use  a  curved  channel  with  an  aspect  ratio  of  6.  Studies  dealing  with 
the  effects  of  heat  transfer  on  turbulent  boundary  layers  include  ones  described  by 
Gibson,  Verriopoulos  and  Nagano  [Ref.  5],  Simon  and  Moffat  [Ref.  6]  and  Chang, 
Humphrey  and  Modavi  [Ref  7], 
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No  studies,  other  than  the  present  one  are  known  to  the  present  authors,  which 
consider  the  effects  of  centrifugal  instabilities  and  vortex  flows  on  local  heat  transfer 
distributions  in  laminar  and  transitional  flow  in  a  channel  with  40  to  1  aspect  ratio. 
Understanding  such  curved  flows  is  important  to  the  design  of  devices  such  as  heat 
exchangers  and  gas  turbine  blades. 

B.  OBJECTIVES 

The  objective  of  this  thesis  is  to  provide  experimental  evidence  as  to  how  curvature 
affects  local  heat  transfer  distributions  on  the  surfaces  of  a  curved  channel  for  Dean 
numbers  from  50  to  200.  This  is  accomplished  by; 

1 .  Qualifying  experimental  procedures,  measurement  techniques  and  channel 
facilities  by  conducting  several  baseline  comparisons  on  the  straight  portion  of  the 
channel. 

2.  Provide  accurate  Nusselt  number  data  for  both  surfaces  of  the  180  degree  curved 
channel,  thus  providing  information  as  to  how  centrifugal  instabilities  and  the  resulting 
vodex  pairs  affect  local  heat  transfer  distributions. 

C.  ORGANIZATION 

Subsequent  to  this  introduction,  Chapter  II  details  the  design  and  construction  of  the 
curved  heated  channel  used  in  this  study.  Chapter  III  discusses  the  experimental 
procedures  and  measurment  techniques.  Chapter  IV  provides  a  series  of  baseline 
qualification  comparisons  and  then  presents  the  results  of  the  Nusselt  number 
measurements.  Chapter  V  presents  conclusions  and  recommendations. 
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II  EXPERIMENTAL  FACILITIES 

A .  CURVED  HEATED  CHANNEL 

The  apparatus  used  in  this  study  is  a  rectangular  curved  channel  shown  in  Figure  1 . 
Within  this  facility,  local  surface  temperatures  and  other  quantities  are  measured  so  that 
the  effects  of  curvature  on  local  Nusselt  number  distributions  can  be  studied.  To  coincide 
with  work  done  on  a  transparent  channel  used  for  flow  visualization,  shown  in  Figure  2, 
identical  internal  dimensions  are  used  for  the  heat  transfer  channel.  The  channel  has  an 
interior  rectangular  cross  section  of  1 .27  cm  by  50.1  cm,  giving  an  aspect  ratio  of  40  to  1 . 
The  heat  transfer  channel  is  different  from  the  transparent  channel  since  it  allows  for 
thermal  expansion  which  results  from  heating.  This  is  accomplished  by  non-rigid  securing 
of  channel  surfaces  and  by  constructing  much  of  the  channel  on  skids  to  allow  for 
longitudinal  expansion.  A  detailed  description  of  the  design  and  construction  of  tha  channel  is 
presented  by  Hughes  [Ref.  8].  A  brief  summary  of  design  and  construction  details  is  also 
given  here. 

1.  Channel  Construction 

A  photograph  of  the  channel  for  heat  transfer  measurements  is  shown  in  Figure  1 . 
Details  of  the  inlet  section  are  shown  in  Figure  3.  The  rectangular  duct  inlet  section  contains 
an  aluminum  honeycomb  and  three  wire  screens  all  placed  normal  to  the  flow  direction. 
These  devices  reduce  spacial  non-uniformities  in  the  flow.  As  the  air  exits  the  25.4  cm  by 
50.8  cm  rectangular  duct  inlet,  it  enters  a  20  to  1  contraction  ratio  nozzle.The  nozzle  is 
constructed  from  the  same  two  continuous  pieces  of  lexan  as  used  for  the  straight  and  curved 
test  section  walls,  thus  eliminating  a  seam  between  the  nozzle  and  the  channel.  Lexan  is  the 
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commercial  name  of  a  polycarbonate  material.  The  shape  of  the  nozzle  is  designed  using  a 
fifth  order  polynomial  with  respect  to  streamwise  distance  to  ensure  that  the  flow  remains 
laminar  and  unseparated  as  it  enters  the  channel. 

After  exiting  the  inlet  nozzle,  air  first  enters  a  straight  section,  244  cm  (  96  inches) 
in  length  with  interior  dimensions  of  1.27  cm  by  50.8  cm  (  .5  inches  by  20  inches  ).  The 
straight  section  allows  hydrodynamically  fully  developed  channel  flow  to  develope  before 
the  flows  enters  the  curved  section.  The  fluid  then  enters  a  180  degree  curved  channel 
section  with  a  convex  surface  radius  of  59.69  cm  (  23.5  inches  )  and  a  concave  surface 
radius  of  60.96  cm  (  24  inches  ).  Upon  exiting  the  curved  section,  air  then  enters  a  second 
straight  section  with  a  length  of  244  cm  (  96  inches  ). 

On  the  convex  surface  of  the  second  straight  section,  a  5.08  cm  (  2.0  inches  )  wide  slot 
is  present  which  allows  insertion  of  a  probe  into  the  channel  to  measure  outlet  mixed  mean 
temperatures.  The  slot  is  located  approximately  57.1  cm  (  22.5  inches  )  from  the  trailing 
edge  of  the  heater  on  the  curved  convex  wall,  or  19.05  cm  (7.5  inches)  downstream  of  the 
end  of  curvature.  A  support  block  was  constructed  to  support  the  area  around  the  slot.  Brass 
inserts  are  placed  just  inside  the  slot  to  allow  insertion  of  two  thin  foam  strips  which 
adhere  to  the  brass  to  provide  an  air  tight  seal  between  the  channel  interior  and  exterior  as 
the  probe  is  traversed.  As  flow  leaves  the  second  straight  portion  it  passes  through  four 
screens,  a  honeycomb,  a  diffuser  with  a  total  angle  of  3  degrees,  and  finally  into  Plenum  #1. 

The  channel  is  designed  to  allow  for  thermal  expansion  in  longitudinal  and  transverse 
directions  which  occurs  as  test  surfaces  are  heated.  As  is  illustrated  in  Figure  7  [Ref.  8], 
the  surfaces  of  the  channel  are  not  rigidly  secured,  but  are  clamped  onto  the  sidewalls  to 
allow  for  such  motion.  In  addition,  both  the  inlet  and  outlet  sections, including  Plenum  #1, 
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are  constucted  on  two  plexiglass  skids  to  allow  for  longitudinal  expansion. 

2.  Blower  Assembly 

The  blower  and  plumbing  connecting  it  to  the  channel  are  shown  schematically  in 
Figure  4  [Ref.  8].  The  device  itself  is  an  ICG  Industries  type  10P  blower,  capable  of 

producing  10.2  cm  of  water  vacuum  at  4.82  m3  /min  volumetric  flowrate.  It  is  used  to 
depressurize  outlet  Plenum  #2  to  pressures  just  below  atmospheric  pressure.  In  doing  so, 
air  is  then  drawn  from  Plenum  #1  through  a  globe  valve  and  a  5.08  cm  inside  diameter  PVC 
pipe.  With  this  system,  flowrates  can  be  varied  to  produce  Dean  numbers  from  10  to  435. 
Flowrates  are  measured  by  means  of  a  pressure  drop  across  a  3.51  cm  (1 .5  inch)  orifice 
plate.  Details  concerning  the  orifice  plate  are  given  in  Hughes  [Ref.  8]  and  also  discussed  in 
Chapter  3. 

To  reduce  vibrations  from  the  blower  to  the  test  section  and  to  minimize  disturbances 
to  the  flow,  a  damper  is  used  to  connect  the  suction  side  of  the  blower  to  Plenum  #2  and  a 
second  damper  is  used  to  connect  Plenum  #2  to  the  orifice  plate.  Both  dampers  consist  of 
thin  sheets  of  plastic  wrapped  around  the  piping  with  air  tight  seals.  The  plastic  is  attached 
at  both  ends  by  double  sided  tape  and  then  tightly  secured  with  standard  hose  clamps  in  order 
to  prevent  any  leaks.  In  addition,  rubber  mounts  are  used  to  mount  the  blower  and  blower 
frame  to  further  reduce  vibrations.  A  photograph  of  this  assembly  is  shown  in  Figure  5. 

3.  Heaters 

Four  etched  foil  heaters,  manufacured  by  the  Electrofilm  Corporation,  are  installed 
on  each  of  the  two  lexan  surfaces,  as  shown  in  Figure  6  [  Ref.  8].  In  the  present  study,  these 
heaters  provide  a  constant  heat  flux  to  both  the  concave  and  convex  surfaces  of  the  channel, 
all  at  the  same  heat  rates,  to  produce  uniform  boundary  conditions  at  the  channel  walls.  The 
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dimensions  of  each  heater  are  38.1  cm  by  152.4  cm  (15.0  inches  by  60.0  inches),  and 
power  capacity  is  2  KW.  Drawings  indicating  how  the  heaters  are  attached  are  shown  in 
Figure  7  and  Figure  8  [Ref.  8]  and  are  discussed  further  by  Schwartz  [Ref.  9j.  The  leading 
edge  of  the  CC1  and  CV1  heaters  are  88.4  cm  from  the  channel  inlet.  The  trailing  edges  of  the 
CC2  and  CV2  heaters  end  approximately  38.1  cm  (15.0  inches)  from  the  completion  of 
channel  curvature.  This  is  about  57.15  cm  (22.5  inches  )  upstream  of  the  location  of  the 
temperature  traverse  probe.  CC1  and  CV1  refer  to  the  two  heaters  on  the  concave  and  convex 
surfaces  of  the  straight  section,  respectively,  and  CC2  and  CV2  refer  to  the  heaters  on  the 
concave  and  convex  surfaces  of  the  curved  section,  respectively. 

Each  heater  is  powered  by  a  Superior  Electric  type  136B  variac  as  shown  in  Figure  9 
[Ref.  8].  With  each  variac,  the  voltage  to  a  heater  may  be  adjusted  between  0  and  115 
volts,  and  the  current  may  be  adjusted  from  0  to  10  amps.  A  detailed  schematic  of  the  heater 
circuitry  is  given  by  Hughes  [Ref.  8]  and  also  shown  in  Figure  1 0  .  Power  inputs  to  each 
heater  are  determined  by  measuring  voltage  drops  across  50  mv,  1 0  amp  shunt  resistors, 
as  well  as  across  each  heater.  The  shunt  resistance  measurement  allows  for  the  heater 
current  to  be  determined.  Details  in  adjusting  desired  heat  input  levels  are  given  in  Chapter 
3. 

4 .  Insulation 

To  minimize  the  loss  of  conducted  heat  to  the  surroundings,  and  maximize  the  amount 
of  heat  into  each  heater  which  is  convected  away,  the  outside  of  the  entire  heated  portion  of 
the  channel  is  insulated  with  black  foam  insulation  manufactured  by  the  Halstead  Company.  A 
detailed  drawing  showing  how  various  layers  of  insulation  are  attached  to  the  channel  to 
produce  a  total  thickness  of  6.0  cm  (2.36  inches)  on  each  side  is  illustrated  in  Figure  7 
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[Ref.  8].  For  this  study,  an  additional  row  of  2.5  cm  (1.0  inch)  insulation  is  attached  along 
all  edges  of  the  heated  channel  to  further  minimize  conduction  from  heated  channel  walls. 
Photographs  of  the  insulation  are  shown  in  Figure  11  and  Figure12. 

5.  Thermocouples 

200  copper-constantan  thermocouples  are  placed  on  channel  surfaces  to  allow 
detailed  spatially  resolved  surface  temperature  measurements  to  be  made.  Figure  8  [Ref.  8] 
shows  the  approximate  locations  of  each  thermocouple.  The  exact  locations  are  listed  in 
Appendix  A.,  as  determined  by  measurements  conducted  as  a  part  of  this  study. 

Measurements  are  based  on  a  coordinate  system  shown  in  Figures  13  and  14  .  Here  X  is  in 
the  streamwise  direction  measured  from  the  inlet  of  the  straight  section,  Z  is  in  the 
spanwise  direction  and  measured  from  the  centerline,  and  Y  is  normal  to  channel  surfaces  in 
the  vertical  direction  measured  from  the  channel  midplane. 

On  the  straight  section, 15  thermocouples  are  placed  on  each  of  the  CC1  and  CV1  surfaces 
in  5  spanwise  rows  of  3  per  row.  The  first  row  is  located  103.64  cm  (  40.8  inches  )  from 
the  inlet  or  15.24  cm  (  6  inches  )  from  the  leading  edge  of  the  heater.  Each  additional  row  is 
spaced  30.48  cm  (  12  inches  )  apart  in  the  X  direction.  The  remaining  170  thermocouples 
are  located  on  the  curved  section  walls  in  5  spanwise  rows  of  17  per  row.  Again,  each  row 
is  spaced  30.48  cm  (12  inches  )  apart  in  the  X  direction.  Of  the  17  thermocouples  in  each 
spanwise  row,  two  are  located  near  the  edges  of  the  channel  with  the  remaining  15  placed 
over  a  spanwise  length  of  approximately  5.08  cm  (  2.0  inches  )  on  each  side  of  the 
centerline.  The  purpose  of  these  thermocouples  is  to  provide  an  accurate  average  local 
Nusselt  number,  as  well  as  a  means  of  determining  local  Nusselt  numbers  in  the  presence  of 
two  pairs  of  Dean  vortices  [Ref.  8). 
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Four  additional  thermocouples  are  used  to  measure  ambient  temperatures  around  the 
outside  of  the  channel,  mixed  mean  temperature  at  the  air  inlet  and  mixed  mean  temperature 
at  the  channel  outlet.  Thermocouple  (TC)  #200  is  placed  just  in  front  of  the  inlet  and 
provides  the  mixed  mean  inlet  temperature.  TC  #201  and  TC  #202  are  placed  respectively 
on  the  outsides  of  the  CC2  and  CV2  surfaces  to  provide  local  ambient  temperatures.  Finally, 
TC  #203  is  used  in  the  temperature  traverse  probe  to  provide  measurement  of  local 
temperatures  used  to  calculate  outlet  mixed  mean  temperatures. 

Each  of  the  204  thermocouples  is  attached  to  a  data  acquisition  channel  and  the 
corresponding  channel  numbers  are  also  listed  in  Appendix  A.  Briefly,  thermocouples  1  to 
30  correspond  to  acquisition  channels  0  to  29,  thermocouples  31  to  115  to  channels  30  to 
79  and  100  to  134,  and  finally,  thermocouples  116  to  204  correspond  to  channels  135  to 
179  and  200  to  243. 

B.  DATA  ACQUISITION 

1  .  Surface  Temperature  Measurement 

Voltages  from  the  204  thermocouples  are  read  by  Hewlett-Packard  type  T20  relay 
multiplexer  card  assemblies  for  T  type  thermocouples.  These  assemblies  are  installed  in  a 
HP3497A  low-speed  Data  Acquisition  /Control  Unit  and  a  HP3498A  Extender.  This  system 
provides  thermocouple  compensation  electronically  such  that  voltages  for  type  T 

thermocouples  are  given  relative  to  0°  C.  This  system  is  connected  to  a  Hewlett-Packard 
9836S  computer  which  processes  voltages  from  the  204  thermocouples  which  are  then 
recorded  into  data  files  along  with  corresponding  temperatures  .  A  schematic  of  this  system 
is  shown  in  Figure  15  [Ref.  8]. 
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2.  Outlet  Mixed  Mean  Temperature  Measurement 

The  acquisition  system  used  for  exit  mixed-mean  temperature  measurments  is 
illustrated  in  Figure  1 6  [Ref.  8}.  The  thermocouple  probe  used  to  measure  the  mixed  mean 
temperature  at  the  exit  of  the  heated  test  section  consists  of  a  thermocouple  attached  to  a 
metal  support.  This  support  is  then  secured  to  a  traversing  block.  An  automated  two- 
dimensional  traversing  device,  which  allows  both  spanwise  and  radial  movement,  is  used  to 
traverse  the  temperature  probe  and  shown  in  Figure  1 7.  The  movements  of  the  motors  on 
the  traversing  devise  are  controled  by  signals  from  the  HP9836S  computer,  which  are  sent 
to  a  Modulynx  Mitas  PMS085-CZAR  Drive  Controller  and  then  to  a  Modulynx  Mitas  PMS085- 
DOJO  Motion  Drive.  The  motor  Drive  sends  signals  to  each  of  the  Superior  Electric  M092- 
FD310  stepping  motors,  which  rotate  drive  screws  which  then  move  two  traversing  blocks. 
The  threads  of  the  drive  screws  are  such  that  one  rotation  gives  a  vertical  or  lateral  probe 
movement  of  .127  cm  (  .05  inches  ).  For  each  traverse,  320  local  temperatures  are 
recorded  over  a  channel  cross  section  area  of  5.08  cm  (  2.0  inches  )  in  the  spanwise 
direction  by  1.02  cm  (  .4  inches)  in  the  radial  direction  in  increments  of  0.127  cm  (  .05 
inches).  Upon  completion  of  the  traverse,  the  controller  is  programmed  to  return  the  probe 
to  its  original  position.  Voltage  readings  from  the  traversing  thermocouple  are  processed  in 
the  same  manner  as  the  channel  surface  thermocouples. 
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III.  EXPERIMENTAL  PROCEDURES 


To  obtain  heat  transfer  distributions,  flow  is  first  induced  in  the  channel,  heater 
power  levels  are  adjusted  to  provide  a  constant  surface  heat  flux  boundary  condition, 
voltages  from  thermocouples  are  read,  surface  temperatures  are  determined,  and  finally 
Nusselt  numbers  are  calculated.  Flow  is  laminar  for  all  tests.  Pairs  of  counter  rotating  Dean 
vortices  are  present  over  a  large  range  of  the  Dean  numbers  studied,  especially  at  farther 
downstream  locations. 

Experimental  procedures  are  discussed  in  three  parts.  In  part  one,  means  to  set 
flowrate  and  Dean  number  are  discussed  along  with  details  on  the  adjustment  of  thermal 
boundary  conditions.  Part  two,  processing  of  surface  temperature  measurements  and  the 
determination  of  local  Nusselt  numbers.  Finally,  part  three  discusses  the  use  of  the 
temperature  traverse  probe,  and  method  used  to  determine  local  and  exit  mixed-mean 
temperatures. 

A.  CHANNEL  PREPARATION 

1.  Flowrate  and  Dean  Number  Setting 

Local  Nusselt  numbers  are  determined  at  Dean  numbers  ranging  from  50  to 
200  in  increments  of  25  using  the  software  program  Nuscurv.  The  first  step  in  each  run  is 

to  set  the  flowrate  to  obtain  the  desired  Dean  number.  The  mass  flowrate,  rh  ,  is 
determined  from  the  equation: 


10 


Where: 


De  =  Dean  Number 
P  =  Air  density 
Acf.  -  Channel  area 

ri  =  Inside  radius  of  channel  curvature 

d  =  Channel  height 

v  =  Kinematic  viscosity 

The  flowrate  is  set  by  regulating  the  pressure  drop  across  a  3.81  cm  (1.5  in )  standard 
AS  ME  orifice  plate,  located  in  piping  between  the  two  exit  plenums.  The  initial  value  of  this 
pressure  drop  is  estimated  using  results  such  as  the  ones  presented  in  Figure  18.  The 
pressure  drop,  in  inches  of  water,  is  then  also  measured  using  a  digital  monometer. 
Afterwards,  an  iteration  process  is  used  to  determine  the  true  Dean  number. 

As  the  iteration  process  is  implemented,  the  pipe  Reynolds  number.  Rep,  is 

calculated  from  the  mass  flowrate,  fn ,  using  the  equation: 

_  rndp 
p  ~  PvAp 


Where: 

dp  =  Pipe  diameter 

Ap  =  Pipe  cross-sectional  area 

A  connected  mass  flowrate  is  then  calculated  using  the  equation: 


m  =  KA^Y^JTpIP 

Where: 

K  =  Flow  coefficient  ( A  function  of  Reynolds  number  which  is  interpolated  from  the 
ASME  Tables  [Ref.  10]) 


Y  =  Expansion  coefficient  (Calculated  from  Holman  and  Gajda  [Ref.  11] 

=  Orifice  Area 

A  P  =  Pressure  drop  across  orifice 

A  corrected  Reynolds  number  is  then  determined  and  the  iterative  process  is  continued  until 
the  mass  flowrate  converges  to  within  one  percent  for  two  successive  calculations.  The 
final  Dean  number  is  then  determined  using  the  equation: 


De  = 


2.  Constant  Heat  Flux  Setting 

The  four  heaters,  CC1,  CC2,  CV1,  and  CV2  are  designed  to  provide  a  constant 
heat  flux  over  both  the  concave,  convex,  and  straight  channel  surfaces.  The  magnitude  of 
the  heat  flux  is  set  based  on  two  contraints.  First,  heat  fluxes  must  be  large  enough  to 
obtain  significant  differences  between  surface  temperatures  and  local  mixed  mean  fluid 
temperatures.  The  second  constraint  limits  maximum  surface  temperatures  for  any  location 
in  the  channel  to  60  Deg.  C.  This  is  necessary  because  the  bonding  agent  used  to  attach 
lexan  sheets  to  cross  beams  is  not  dependable  at  temperatures  above  65  Deg.  C.  For  most 
runs,  the  maximum  temperature  occurs  at  TC  #199,  which  is  located  in  the  final  spanwise 
row  of  the  convex  wall  of  the  curved  section.  A  curve  showing  heat  flux  required  per  heater 
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to  reach  a  maximum  temperature  of  60  Deg.  C,  obtained  from  experimental  tests,  is  shown 
in  Figure  19  for  Dean  numbers  from  50  to  200. 

To  determine  the  heat  flux,  the  voltage  drops  across  each  shunt  and  each  heater  are 
measured  using  a  simple  volt-ohm  meter  and  are  entered  into  a  program  called  Temcurv. 
Figure  20  shows  voltage  drops  across  shunt  resistors  which  are  needed  to  produce  desired 
power  levels  in  each  heater.  Equations  relating  these  voltage  drops  to  power  are  given  by: 


and 


1=  AVS/R 

P=  A  VH  x  I 


Where: 

I  =  Current  (amps) 

P  =  Power  (watts) 

A  Vs  =  Voltage  drop  across  the  shunt  (volts) 

A  Vh  =  Voltage  drop  across  the  heater  (volts) 

R  =  Shunt  resistance  (ohms) 

B.  SURFACE  HEAT  TRANSFER  DISTRIBUTION  MEASUREMENTS 
1.  Surface  Temperature  Measurements 

The  system  is  heated  for  approximately  8-9  hours  prior  to  data  acquisition  to 
reach  steady  state  and  to  insure  that  thermal  equilibrium  exists.  Great  care  is  taken  to  insure 
that  maximum  surface  temperatures  never  exeed  60  Deg.  C.  by  routinely  checking 
temperatures  using  the  temperature  acquisition  program,  Tcheck.  Figure  21  provides  a 
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curve  from  Run  #120689.2125  which  illustrates  how  temperature  tends  to  increase  with 
time.  When  steady  state  conditions  are  reached,  voltage  readings  from  each  thermocouple 
are  recorded  using  Temcurv  and  converted  into  Deg.  C  using  calibration  equations 
determined  by  Hughes  [Ref.  8].  For  thermocouples  1  through  200,  201, 203  and  204, 
voltages  are  converted  into  Deg.  C  using  the  equation: 

T  =  -1.988321  +  29436.53V  -  2622880  V2  +  36932000 V3 

For  thermocouple  202  the  equation  is: 

T  =  -1.816347  +  2928.91V  -  2622880V2  +  373752000 V3 

Where: 

T  =  Temperature  Deg.  C 
V  =  Thermocouple  voltage  (volts) 

These  equations  are  utilized  in  the  software  program  Nuscurv3B. 

The  surface  thermocouples  are  attached  behind  Lexan  surfaces  exposed  to  the 
flow  field.  Consequently,  there  is  a  thermal  contact  resistance  between  thermocouples  and 
convective  Lexan  surfaces.  The  temperature  difference,  A  Tw^,  between  the  actual 
surface  and  that  measured  by  the  thermocouple  is  given  by  the  equation: 

A  Twcon  =  q  conv  (CR) 

Where: 

q'  conv  =  Convected  heat  flux  per  heater  area 
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CR  =  Thermal  contact  resistance  (°Cm2/ watt) 

Experimental  procedures  used  to  determine  CR  are  described  by  Bella  [Ref.  12].  From  his 
work,  CR  equals  3.4  x  10-3  o  cm2  /  watt. 

2.  Nusselt  Number  Calculations 

The  program  Nuscurv3B  is  used  to  determine  local  Nusselt  numbers.  The 
process  begins  with  determination  of  the  convective  power  transfer  using: 

Q  conv  —  Q  in  "  Q  cond  ~  Q  ad 


Where: 

q  conv  =  Power  convected  to  the  channel  air  flow 
q  m  =  Power  to  the  heaters  (V  x  I) 
q  cond  =  Power  loss  by  conduction 
Qb i  =  Power  loss  by  radiation,  set  equal  to  zero 

Power  loss  by  conduction  is  given  by  Bella  [Ref.  12],  who  conducted  tests  on  a 
full  scale  prototype  model  of  the  curved  heated  section.  His  results  show  that  conduction 
losses  are  given  by: 

Q  cond  =  -0.0294  +  0.4222  (Twcorr  "  Tambienl)  *  0.0015  (Twcorr  -  Tajnbient)^ 


Where: 

Twcorr  =  Corrected  local  surface  temperature 
Tambicnt  =  Ambient  temperature  for  surroundings 


The  convective  heat  flux  is  then  determined  using: 


Q  conv 


=  Q  conv /A 


sec 


Where: 

Asa;  =  Surface  area  for  each  of  the  four  heater  surfaces, CC1,CC2,CV1  and  CV2 

With  the  convective  heat  flux  known,  the  local  mixed  mean  temperature  at  any 
streamwise  channel  location  is  determined  using: 

hnm(x)  =  hnmin  conv(x)  b  X  ]  /  Cp  (  ifl  ) 


Where: 

Wn(x) 

tirunin 

b 

x 

Q  "conv 

Cp 


Local  mixed  mean  Temperature 
Inlet  mixed  mean  temperature 
Spanwise  width 
Streamwise  distance 

Heat  flux  over  streamwise  distance  x 
Specific  heat  of  air 


Local  heat  transfer  coefficients  are  then  determined  using  the  equation: 


h  —  Q  conv  /  (  Twcorr  -  tmm  ) 
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Finally,  the  local  Nusselt  number,  Nu,  is  calculated  from: 


Nu  =  h  ( DH)  /  k 


Where: 

Dh  -  Hydraulic  diameter  ( 2  x  channel  width  =  .0254  m) 
k  =  Thermal  conductivity  of  air  at  200  Celcius  ( .02563  W  /  m  K  ) 

C.  EXIT  MIXED  MEAN  TEMPERATURE  MEASUREMENTS  AND 
DETERMINATION 

1 .  EXIT  Mixed-Mean  Temperature  Determination 

A  two  step  procedure  is  used  to  determine  the  local  mixed-mean  temperature 
from  the  convective  heat  flux.  This  temperature  is  determined  at  the  location  of  the  two- 
dimensional  traverse  slot,  which  is  57.1  cm  from  the  trailing  edge  of  the  heater  on  the 
convexly  curved  wall.  First,  the  mixed-mean  temperature  is  estimated  at  the  downstream 
edge  of  the  heaters  using  procedures  described  in  the  previous  section.  This  temperature  is 
noted  as,  tm^, .  In  the  second  step,  an  energy  balance  is  then  used  to  estimate  the  heat  loss 
from  the  flow  to  the  unheated  portion  of  the  channel,  which  extends  from  the  trailing  edge 
of  the  heater  to  the  traversing  slot.  This  energy  balance  is  given  by: 

d?  Cptnuru  —  ITI  Cpvnirno  +  Atoth[(tnimi  tmmo)/2  -  T^v] 

Where: 

tmrr .  =  Local  mixed-mean  temperature  at  the  traverse  slot. 

Tw  =  Temperature  of  the  unheated  channel  surface  (approximately  equal  to  ambient 
temperature). 


Atot  =  Concave  and  convex  surface  area  for  this  portion  of  the  channel, 
h  =  Heat  transfer  coefficient  at  last  span  wise  row  of  thermocouples 
This  equation  is  employed  in  the  software  plotting  program  Tmxdmean,  in  the  form 
where  slot  temperature  t^no  is  given  by: 


*mmo 


KAJ 


mCp  + 


hAtot 

2 


2 .  Exit  Mixed-Mean  Temperature  Measurement 

After  the  channel  is  thermally  steady  state,  TC  #204  is  set  in  place  into  the 
temperature  traversing  carriage.  Using  the  program  Temtrav,  320  voltage  readings  are 
taken  over  a  5.08  cm  (2.0  inch)  by  1.02  cm  (0.4  inch)  area  as  the  probe  moves  in  the 
traverse  slot.  Voltage  readings  are  transformed  into  Deg.  C  using  the  equation  described 
the  previous  section.  The  mixed  mean  temperature  is  then  given  by: 


=  Aj7$U,dA 


Where: 


U  = 


Mixed  mean  temperature 
Channel  traversed  cross  section  area 


u 


Bulk  velocity  (spacially  averaged) 
Local  mean  stream  wise  velocity 
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t  =  Local  mean  temperature 

The  software  program  Tmxdmean  is  used  to  determine  the  mixed-mean  temperature.  This 
is  accomplished  using  an  approach  which  assumes  that  the  area  corresponding  to  each 
temperature  measurement  is  the  same: 


1 

"  320 


Velocity  data,  U  and  u; ,  were  obtained  using  a  miniature  five-hole  pressure  probe 
described  by  Baun  [Ref.  13].  The  velocity  data  used  to  determine  mixed-mean  temperatures 
in  the  present  study  are  the  same  as  given  by  Hughes  [Ref.  8]  and  Tuzzolo  [Ref.  14]. 
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IV.  RESULTS  AND  DISCUSSION 

The  objective  of  the  present  channel  study  is  to  examine  the  effects  of  curvature  and  the 
resulting  centrifugal  instabilities  on  local  heat  transfer  distributions.  Rows  at  Dean  numbers 
ranging  from  50  to  200  are  investigated.  The  results  of  these  tests  are  discussed  in  four 
parts.  First,  two  types  of  baseline  tests  are  discussed,  which  qualify  measuring  procedures 
and  apparatus.  This  section  includes  discussion  of  spanwise-averaged  Nusselt  numbers 
from  the  straight  portion  of  the  channel  as  well  as  comparison  of  measured  exit  mixed- 
mean  temperatures  with  values  determined  from  energy  balances.  Second,  spanwise- 
averaged  Nusselt  number  distributions  from  concave  and  convex  surfaces  are  presented. 
Third,  the  effects  of  counter-rotating  Dean  vortex  pairs  on  the  local  Nusselt  number 
distributions  are  detailed.  Finally,  in  the  fourth  section,  surveys  of  local  temperatures  used 
to  determine  the  mixed  mean  temperature  are  presented  in  conjunction  with  mean  velocity 
data  obtained  by  Hughes  [Ref.  8]  and  Tozzolo  [Ref.  14]. 

A.  BASELINE  RESULTS 

1.  Spanwise-Averaged  Nusselt  Numbers  From  the  Straight  Portion  of 
the  Channel. 

The  walls  of  the  straight  portion  of  the  channel  are  heated  to  insure  that  flow  is 
thermally  and  hydrodynamically  fully  developed  prior  to  entering  the  curved 
section.Because  this  portion  of  the  channel  is  straight,  much  information  is  available  on  the 
behavior  and  distributions  of  local  Nusselt  numbers.  First,  because  the  flow  is  symmetic  in 
the  direction  normal  to  the  walls.  Nusselt  numbers  should  be  the  same  on  the  top  and 
bottom  walls  at  each  streamwise  location.  Second,  the  analytical  solution  for  flow  between 
two  infinite  parallel  plates  with  constant  heat  flux  boundary  conditions  gives  the  fully 
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developed  Nusselt  number  value.  Third,  a  numerical  solution  is  used  to  predict  thermal 
entry  length  behavior.  All  three  of  these  are  now  discussed  in  regard  to  measured  Nusselt 
number  distributions. 

a.  Comparison  of  Results  from  the  Straight  CC1  and  CV1  Surfaces 
Figure  22  and  Figure  23  show  span  wise-averaged  Nusselt  numbers  for  Dean 

numbers  of  100  and  125  for  the  CC1  and  CV1  surfaces,  respectively.  The  graphs  show 
that  the  Nusselt  number  distributions  are  nearly  identical  over  the  entire  lengths  of  the 
straight  sections,  from  x/d  of  82.1  cm,  which  is  the  location  of  the  first  spanwise  row  of 
thermocouples,  to  x/d  of  178.6  cm,  the  beginning  of  curvature. 

b.  Comparison  of  Experimental  Data  with  the  Analytic  Solution  for 
Thermally  Fully  Developed  Flow 

Kays  and  Crawford  [Ref.  15],  indicate  that  the  average  Nusselt  number  is  8.24 
for  fully  developed  laminar  flow  between  two  infinite  parallel  plates  with  constant  heat  flux 
boundary  conditions.  Figures  22  and  23  show  measured  Nusselt  numbers  ranging  from 
8.45  to  8.1 1,  which  agree  excellently  with  the  value  of  8.24,  for  x/d  from  130.4  to  178.6. 
As  indicated  earlier,  these  data  are  obtained  for  Dean  numbers  of  100  and  125  and  are 
approximately  constant  with  x/d  only  after  thermal  entry  length  effects  disappear  and  the 
flow  becomes  thermally  fully  developed. 

c.  Comparison  of  Thermal  Development  Effects  with  Numerical 
Solutions 

Span  wise-averaged  Nusselt  numbers  for  Dean  Numbers  of  100  and  125  are 
compared  to  results  from  numerical  predictions  in  Figures  24  and  25.The  numerical 
prediction  results  are  obtained  using  the  STAN5  boundary  layer  program,  described  by 
Kays  and  Crawford  [Ref.  15].  Exact  experimental  conditions  are  used  as  initial  conditions 
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and  boundary  conditions,  and  include  the  effects  of  the  unheated  starting  length  between 
the  nozzle  exit  and  heater  leading  edges.  The  data  and  predictions  show  excellent  agreement 
along  the  entire  length  of  the  straight  portion  of  the  channel,  from  x/d  of  70.0  to  178.64, 
which  validates  measurements  made  in  the  thermal  entry  length  region  of  the  channel. 

Predicted  Nusselt  numbers  for  the  straight  portion  of  the  channel  for  all  Dean 
Numbers  studied  are  given  in  Figure  26.  At  each  Dean  number,  Nusselt  numbers  are 
infinite  at  the  initial  heating  point  With  stream  wise  development,  Nusselt  numbers  then 
decrease.  At  each  particular  streamwise  location,  spanwise-averaged  Nusselt  numbers 
increase  with  Dean  number  because  the  entry  length  needed  for  flow  to  become  thermally 
fully  developed  increases.  The  same  Nusselt  numbers  are  non-dimensionalized  in  Figures 
27  and  28.  Predicted  and  measured  non-dimensional  Nusselt  numbers  for  the  concave 
surfaces  are  compared  in  Figures  29  and  30.  Figures  31  and  32  then  give  the  same  type  of 
infomation  for  the  convex  surface. 

2.  Mixed-Mean  Temperature  Comparison 

Figures  33  through  39  present  mixed-mean  temperature  distributions  along  the 
length  of  the  channel,  which  are  determined  from  the  convected  heat  flux  into  the  flow. 
These  figures  also  compare  these  distributions  to  mixed-mean  temperatures  determined 
from  local  measurements  of  temperature  and  velocity  obtained  from  probes  placed  in  a  slot 
located  downstream  of  the  last  heater.  In  figures  33-39,  results  are  given  for  Dean 
numbers  of  58.5,  73.7, 98.1,  126.0,  149.8,  173.5  and  198.4.  In  all  cases,  mixed-mean 
temperatures  measured  at  the  exit  of  the  curved  portion  closely  match  values  at  x/d  of  356.6 
from  the  energy  balance.  This  further  validates  experimental  procedures,  experimental 
apparatus,  and  the  energy  balances  employed,  and  confirms  that  the  conduction  loss 
estimates  used  in  determining  the  convected  heat  fluxes  are  accurate.  Small  differences  in 
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the  two  exit  mixed-mean  temperatures  sometimes  result  from  small  uncertainties  in 
measured  temperatures  and  velocities.  Mean  temperature  distributions,  mean  velocity 
distributions  and  distributions  of  the  mean  temperautre/mean  velocity  product  used  to 
determine  the  local  mixed-mean  temperature  are  presented  in  Part  D.  of  Chapter  4. 

3.  Baseline  Data  Summary 

All  baseline  test  results  validate  experimental  procedures,  the  measurement 
approach,  the  energy  balances  employed,  and  the  facilities  used  in  this  study,  and  provide 
confidence  that  the  measurements  taken  in  the  curved  portion  of  the  channel,  located  just 
downstream  of  the  straight  portion,  represent  actual  flow  behavior. 

B .  SPANWISE-AVERAGED  NUSSELT  NUMBER  DISTRIBUTIONS 
In  this  section,  spanwise-averaged  Nusselt  numbers  are  presented,  which  are  non- 
dimensionalized  using  the  spanwise-averaged  Nusselt  number  at  x/d  of  130.4.  This  is  done 
to  minimize  some  of  the  scatter  in  dimensional  Nusselt  numbers  measured  at  Dean  numbers 
lower  than  100  and  higher  than  125. 

Figures  40, 41, 42  and  43  present  Nusselt  number  ratios  for  concave  and  convex 
surfaces  at  seven  Dean  Numbers.  At  x/d  values  between  0  and  82.1,  no  Nusselt  numbers 
are  plotted  because  this  portion  of  the  channel  is  unheated.  Between  x/d  of  82.1  and 
approximately  1 10,  Nusselt  numbers  increase  with  Dean  number,  with  trends  consistant 
with  predictions  which  account  for  thermal  entry  length  effects.  At  x/d  of  about  1 10  to  130, 
flow  is  hydrodynamically  and  thermally  fully  developed  and  Nusselt  numbers  are 
approximately  constant  with  stream  wise  distance.  As  shown  in  Figure  22,  the  spanwise- 
averaged  Nusselt  number  is  about  8.24  for  Dean  numbers  of  100  and  125,  which  is  in 
agreement  with  the  theory  for  the  same  flow.  The  same  behavior  exists  for  Dean  numbers 
up  to  about  150,  whenever  the  flow  is  fully  laminar  and  no  transition  phenomenon  are 
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present.  As  x/d  approaches  160-180,  span  wise-averaged  Nusselt  number  ratios  increase 
for  Dean  numbers  of  175  and  200,  as  a  result  of  initiation  of  transition  phenomena.  Peak 
values  exist  at  x/d  of  about  178.6,  which  is  just  at  the  begining  of  curvature.  Here,  some 
mixing  of  flow  is  present  that  does  not  occur  otherwise.  This  trend  occurs  on  both  concave 
and  convex  surfaces,  but  it  is  more  pronounced  on  the  convex  surface. 

As  x/d  becomes  greater  than  180,  effects  of  curvature  are  evident  in  measurements 
from  both  the  concave  and  convex  surfaces.  These  effects  are  complex  and  vary  not  only 
with  the  Dean  number,  but  also  with  the  surface  considered.  On  the  concave  surface  at 
Dean  numbers  of  50, 75, 100  and  125,  Nusselt  number  ratio  distributions  are  the  same  for 
x/d  of  178.6  to  approximately  226.9.  For  each  case,  the  Nusselt  numbers  decrease 
drastically  as  a  result  of  the  stabilizing  influence  of  convex  curvature,  which  causes  a 
stabilizing  effect  initially  near  the  convex  surface  before  propagating  across  the  width  of  the 
channel. 

At  x/d  greater  then  226.9  for  both  surfaces,  the  presence  of  Dean  vortices  is  evident 
in  local  Nusselt  number  ratios,  which  gradually  increase  with  streamwise  development. 
These  vortices  develop  because  of  the  imposition  of  centrifugal  instabilities.  For  x/d  from 
226.9  to  299.4,  the  lowest  Nusselt  number  ratios  correspond  to  a  Dean  numbers  of  50,  and 
the  highest  correspond  to  Dean  numbers  of  100  and  125.  At  these  higher  Dean  numbers, 
centrifugal  instabilities  and  the  accompanying  Dean  vorices  are  more  intense,  and  thus, 
Nusselt  numbers  are  also  higher. 

At  Dean  numbers  of  150, 175  and  200,  results  from  both  the  concave  and  convex 
surfaces  show  evidence  of  transition  phenomena  resulting  in  the  onset  of  some  chaotic 
motions.  This  is  particularly  so  for  x/d  ranging  from  154.5  to  178.6,  just  prior  to 
curvature,where  Nusselt  number  ratios  are  locally  higher.  At  x/d  from  178.6  to  300,  ratios 
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from  the  convex  surface  are  generally  lower  than  from  the  concave  surface.  This  occurs 
because  of  the  stabilizing  and  destabilizing  influences  of  the  convex  and  concave  surfaces, 
respectively.  Row  visualization  results  for  the  same  experimental  conditions  show 
considerable  mixing. 

C .  EFFECTS  OF  COUNTER-ROTATING  DEAN  VORTEX  PAIRS  ON 

LOCAL  NUSSELT  NUMBER  DISTRIBUTIONS 

Local  Nusselt  numbers  versus  non-dimensional  spanwise  length,  z/d,  for  both  concave 
and  convex  surfaces  at  Dean  numbers  of  58.5,  73.7,  98.1,  126.0,  149.8,  173.5  and 
198.4,  are  given  in  Figures  44-57.  These  data,  especially  those  from  the  concave  surface, 
illustrate  the  effects  of  pairs  of  counter-rotating  Dean  vortices  on  local  Nusselt  number 
distributions. 

In  general,  local  Nusselt  numbers  on  the  convex  surface  are  generally  spanwise 
uniform,  regardless  of  Dean  number.  Exceptions  occur  at  Dean  numbers  of  149.8, 173.5 
and  198.4,  where  some  slight  variations  are  present  probably  due  to  vortex  pairs.  On  the 
concave  surface,  local  Nusselt  numbers  develop  significant  variations  in  the  spanwise 
direction  as  the  flows  develop  downstream  as  a  result  of  vortex  upwash  and  downwash 
regions.  Local  Nusselt  number  values  are  augmented  near  downwash  regions  and 
diminished  in  upwash  regions.  Such  spanwise  variations  coincide  with  increases  of 
spanwise-averaged  Nusselt  numbers,  and  are  most  significant  for  Dean  numbers  of  98.6 
and  greater.  Such  behavior  is  also  consistant  with  flow  visualization  results  which  show 
counter-rotating  vortex  pairs  for  the  same  experimental  conditions. 

Focussing  on  Nusselt  number  results  for  a  Dean  number  of  98.6  and  comparing  them 
to  the  spanwise  averaged  Nusselt  number  data  presented  in  section  B,  interesting 
conclusions  are  made.  Data  at  x/d  of  202.8  to  226.9  show  that  local  Nusselt  numbers  are 
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span  wise  uniform.  These  locations  coincide  with  regions  where  span  wise  averaged  Nusselt 
numbers  are  lower  than  values  at  other  x/d,  as  a  result  of  the  imposition  of  centrifugal 
instabilities.  As  x/d  increases  to  251.1, 275.2  and  299.4,  spanwise-averaged  Nusselt 
numbers  increase  and  local  Nusselt  numbers  show  considerable  span  wise  variations  due  to 
the  formation  of  Dean  vortex  pairs. 

Such  trends  are  present  at  all  Dean  numbers  greater  than  or  equal  to  98.6.  However,  at 
Dean  numbers  of  173.5  and  198.4,  local  Nusselt  numbers  show  spanwise  variations  at  x/d 
as  low  as  202.8  to  226.9. 

D .  LOCAL  MIXED-MEAN  TEMPERATURE  SURVEY  DATA 

1 .  Local  Mean  Temperature  Distributions 

Mean  temperature  surveys  are  presented  for  Dean  numbers  of  58.5, 73.7, 98.1, 
126.0,  149.8,  173.5  and  198.4  in  Figures  58-64.  Here,  mean  refers  to  a  time-average.  For 
Dean  numbers  of  58.5  and  73.7,  spanwise  periodicities  result  due  to  vortex  pairs.  At  Dean 
numbers  from  98.1  to  198.4,  the  local  mean  temperature  distributions  are  fairly  spanwise 
uniform,  with  only  slight  variations  due  to  vortex  pairs.  Such  variations  are  present  at 
lower  Dean  numbers  because  Dean  vortices  form  only  after  a  considerable  amount  of 
streamwise  development  Such  vortices  continue  to  exist  at  higher  Dean  numbers,  but  they 
are  unsteady  over  spanwise  distances  greater  than  the  spacing  between  vortex  pairs,  which 
makes  them  appear  spanwise  uniform  in  the  time-averaged  flow  field. 

2.  Local  Mean  Velocity  Distributions 

Mean  velocity  distributions  are  given  by  Hughes  [Ref.  8]  and  Tuzzolo  [Ref.  14]. 
Results  for  Dean  numbers  of  50, 75, 100  and  150  are  shown  in  Figures  65-68.  Qualitative 
trends  of  the  velocity  contours  coincide  nicely  with  local  mean  temperature  contours; 
velocity  profiles  at  low  Dean  numbers  show  some  spanwise  variations,  while  at  higher 
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Dean  numbers  the  profiles  are  fairly  span  wise  uniform. 

3 .  Distributions  of  the  Mean  Temperature-Mean  Velocity  Product 

Mean  temperature-mean  velocity  product  distributions  are  presented  in  Figures  69- 
75  for  Dean  numbers  of  58.5,  73.7,  98.1,  126.0,  149.8,  173.5  and  198.4.  The  graphs 
show  good  span  wise  uniformity  for  all  Dean  numbers,  except  for  some  of  the  lower 
values,  as  expected.  As  mentioned  in  Part  C  of  Chapter  3,  these  temperature- velocity 
product  data  are  used  to  compute  the  local  mixed  mean  temperatures  discussed  in  Chapter 
4,  Part  A.,  Section  2. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  effects  of  curvature  and  the  resulting  centrifugal  instabilities  on  local  heat  transfer 
distributions  are  studied  in  a  curved  channel  at  Dean  numbers  ranging  from  50  to  200.  The 
channel  is  designed  with  a  rectangular  cross  section  of  1.27  cm  by  50.1  cm  giving  an 
aspect  ratio  of  40  to  1.  These  interior  dimensions  coincide  with  those  of  another  channel 
which  is  transparent  and  used  for  flow  visualization.  In  the  channel  used  for  thermal 
studies,  the  flow  is  heated  in  a  2.44  m  straight  portion,  located  just  prior  to  the  curved 
portion  to  obtain  flow  which  is  hydrodynamically  and  thermally  developed.  The  flow  then 
enters  a  180  degree  heated  curved  section  with  a  convex  surface  radius  of  59.69  cm  and  a 
concave  surface  radius  of  60.96  cm.  Surface  temperatures  and  ambient  temperatures  are 
measured  using  204  thermocouples,  whose  outputs  are  used  to  determine  local  Nusselt 
numbers  as  constant  heat  flux  boundary  conditions  are  maintained  on  all  heated  surfaces. 

Several  baseline  comparisons  qualify  flow  behavior  as  well  as  the  experimental 
procedures  and  measurement  techniques.  In  the  first  of  these,  spanwise-averaged  Nusselt 
numbers  from  the  two  walls  of  the  straight  portion  of  the  channel  match  with  each  other 
very  well.  As  thermal  entry  length  effects  disappear,  dimensional  Nusselt  numbers 
approach  8.24  as  expected  from  the  theory  for  thermally  fully  developed  flow  between  two 
parallel  infinite  plates  with  constant  heat  flux  boundary  conditions.  In  the  thermally 
developing  region,  measured  Nusselt  numbers  match  results  from  numerical  predictions. 
Finally,  measured  mixed-mean  temperatures  at  the  exit  of  the  curved  section  match  values 
determined  from  energy  balances.  These  baseline  comparisons  thus  validate  experimental 
procedures,  the  measurement  approach  and  the  facilities  used,  and  provide  confidence  of 


28 


the  accuracy  of  measurements  from  the  curved  portion  of  the  channel. 

Spanwise-averaged  Nusselt  numbers  are  non-dimensionalized  with  respect  to  the 
Nusselt  number  determined  at  x/d  of  130.4.  For  Dean  numbers  of  50, 75, 100  and  125, 
these  ratios  approach  constant  values  on  the  straight  portion  of  the  channel  for  x/d  from  1 10 
to  180.  At  smaller  x/d,  Nusselt  number  ratios  increase  with  Dean  number  as  a  result  of 
thermal  entry  length  effects.  For  Dean  Numbers  of  150, 175  and  200  increases  in  Nusselt 
number  ratios  are  present  for  160  <  x/d  <  180  as  a  result  of  transition  phenomena. 

The  stabilizing  influences  of  convex  curvature  are  evidenced  at  the  onset  of  curvature 
by  drastic  decreases  in  Nusselt  number  ratios  for  both  surfaces  at  Dean  numbers  of  50, 75, 
100  and  125.  These  effects  initially  produce  a  stablilizing  effect,  which  occurs  primarily 
near  the  convex  surface  and  then  propagates  across  the  width  of  the  channel.  As  x/d 
increases,  the  formation  of  Dean  vortices  as  a  result  of  centrifugal  instabilities  is  evidenced 
by  gradual  increases  of  the  local  Nusselt  number  ratios.  Here,  Nusselt  number  ratios  are 
generally  lower  on  the  concave  surface,  where  the  lowest  Nusselt  number  ratios  are  present 
at  a  Dean  number  of  50  and  the  highest  values  exist  at  Dean  Numbers  of  100  and  125. 
Higher  Nusselt  numbers  are  expected  at  higher  Dean  numbers  because  centrifugal 
instabilities  and  the  resulting  Dean  vortices  are  more  intense.  Higher  Nusselt  numbers  are 
expected  near  the  concave  surface  because  of  the  intensity  of  secondary  flows  from  counter 
rotating  vortex  pairs,  which  are  much  less  intense  near  convex  surfaces. 

At  Dean  numbers  of  150, 175  and  200,  effects  of  transition  phenomena  are 
considerably  more  apparent.  Nusselt  number  ratios  show  local  increases  on  the  straight 
portion  of  the  channel  just  prior  to  curvature.  After  curvature  is  imposed,  trends  are 
generally  similar  ones  observed  at  the  lower  Dean  numbers,  except  that  Nusselt  number 
ratios  are  higher  probably  because  of  increased  flow  mixing  associated  with  transition 


phenomena. 

Local  Nusselt  numbers  from  the  convex  surface  are  fairly  spanwise  uniform, 
regardless  of  Dean  number.  On  the  concave  surface,  however,  significant  spanwise 
variations  are  evident,  especially  at  Dean  numbers  of  98.6  or  greater,  which  result  from 
upwash  and  downwash  regions  in  vortex  pairs.  Such  spanwise  variations  correspond  to 
increases  in  span  wise-averaged  Nusselt  numbers  which  occur  with  stream  wise 
development. 

Mean  temperature  and  mean  velocity  survey  data  show  some  spanwise  periodicity  at  the 
lower  Dean  numbers  due  to  upwash  and  downwash  regions  formed  by  vortex  pairs.  As  the 
Dean  number  increases,  these  surveys  show  distributions  which  are  fairly  spanwise 
uniform. 

B  .  RECOMMENDATIONS 

Recommendations  are,  first,  repeat  the  present  study  at  Dean  numbers  from  200  to  450 
to  obtain  additional  information  on  the  connections  between  centrifugal  instabilities.  Dean 
vortices  and  transition  from  laminar  flow  to  turbulent  flow.  Second,  Nusselt  numbers 
should  be  measured  in  situations  where  the  heat  flux  differs  on  the  two  surfaces  of  the 
channel. 
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APPENDIX  A 

THERMOCOUPLE  PLACEMENT  TABLES 
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APPENDIX  B 
FIGURES 


The  following  pages  contain  the  figures  used  in  the  development  of 
this  thesis. 
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Curved  Heat  Transfer  Channel 
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Curved  Transparent  Channel 


Valve 


Schematic  of  Blower  Assembly 
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Figure  5.  Photograph  of  Blower  Assembly 
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Figure  9.  Channel  Heater  Control  Variacs 


45 


*-  <v  r>  _  * 
-  n  it  n  »»  i*  n 

>i  »•  n  n 


*r  «  •=  • 


OOOOOOOO 

OOOOOOOO 

OOOOOOOO 


C7,C7,C7'0*C7,0'C?0 

nnitnnnnn 


O 


OOOOOOOO 

>>>>>>>> 


O  *-  oj  n  tn  to  r-  a? 


WVW-1 


-VVVVWfc 


V\/WV' 


Nj' 

rWVvV 


I 


V\AAA/ — 

r> 

I 


ft 


5 


f — \AA/VV—9JJ  r — AAAA/'Oj 

r  I  ?  .  I 


_ L 


1 


Figure  10.  Heater  and  Control  Circuit  Wiring  Diagram 
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Figure  1 1 .  Photograph  of  Channel  Insulation 
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Figure  15.  Data  Acquisition  System  For  Surface  Temperature 
Measurement 
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Figure  1 6.  Data  Acquisition  System  For  Mixed-Mean 
Temperature  Measurement 
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PRESSURE  DROP  VS  DEAN  NUMBER 
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Figure  18.  Pressure  Drop  Across  Orifice  Plate  vs.  Dean 
Number 
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Figure  19.  Dean  Number  Versus  Power  Input  Per  Heater 
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VOLTAGE  DROP  VS  POWER 
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Figure  20.  Voltage  Drop  Across  Shunt  Versus  Heat  Flus  Per 
Heater 
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Figure  21.  Temperature  Increase  vs.  Time  For  Run#120689 
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Figure  22.  Spanwise  Average  Nusselt  Number  Distribution  vs. 
x/d  ,  Concave  Surface,  De  =  100  and  125 
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Figure  24.  Spanwise  Averaged  Nusselt  Number  vs.  x/d, 

Comparison  Of  Experimental  And  Analytical  Values 
Concave  Surface,  De  =  100  and  125 
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Figure  25.  Spanwise  Averged  Nusselt  Number  xs.  x/d, 

Experimental  and  Analytical  Values,  Convex  Surface, 
De  =  100  and  125 
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Figure  26.  Predicted  Straight  Channel  Flow  For  Averaged 
Nusselt  Number  vs.  x/d 
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Figure  27.  Predicted  Straight  Channel  Flow  For  Non- 
Dimensional  Nussett  Number  vs.  x/d 
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Figure  28.  Predicted  Staight  Channel  Flow  For  Non-Dimensional 
Nusselt  Number  vs.  x/d 
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Figure  30.  Nu/Nuo  vs.  x/d,  Concave  Surface,  Experimental  and 
Predicted  Values  For  De  =  150,  175  and  200 
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Figure  33.  Mixed-Mean  Temperature  vs.  x/d  Comparison, 
De  =  58.5 
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Figure  40.  Nu/Nuo  vs.  x/d,  Concave  Surface,  All  Dean  Numbers 
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Figure  47.  Spanwise  Nusselt  Number  Distributions  vs  z/d  For 
Curved  Convex  Surface,  De  =  74.2 
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Figure  50.  Spanwise  Nusselt  Number  Distributions  vs.  z/d  For 
Curved  Concave  Surface,  De  =  126.4 
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Fiqure  59.  Measured  Exit  Local  Mean  Temperature  Contours. 
De  =  73.7 
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Figure  60.  Measured  Exit  Local  Mean  Temperature  Contours. 
De  =  98.1 
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Figure  61 .  Measured  Exit  Local  Mean  Temperature  Contours. 
De  =  126.0 
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Figure  63.  Measured  Exit  Local  Mean  Temperature  Contours. 
De  =  173.5 
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Figure  64.  Measured  Exit  Local  Mean  Temperature  Contours. 
De  =  198.4 
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Figure  65.  Mean  Velocity  Contours,  De  =  50 
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Figure  66.  Mean  Velocity  Contours,  De  =  75 
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Figure  68.  Mean  Velocity  Contours,  De  =  150 
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Figure  69.  Mean  Temperature/Mean  Velocity  Product  Contours, 
De  =  58.5 
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Figure  70.  Mean  Temperature/Mean  Velocity  Product  Contours, 
De  =  73.7 
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Figure  71 .  Mean  Temperature/Mean  Velocity  Product  Contours, 
De  =  98.1 
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Figure  72.  Mean  Temperature/Mean  Velocity  Product  Contours, 


CO 

* 

CD 

\r 


0 

P 


co 

CD 


* 

(3 

cn 

■r— 4 

C\J 

*--4 


z> 

CK 


i 

CD  CD  ^  CM  S 


p/A 


in 

ui 

CD 


rvj 

co 


cm  r» 
cn  CS3 


cr 

OC 


1A 

\ 

£ 


oj  rs.  oj  OJ 
ro  id  fs.  cn 


in  o?  cn 
c n  pi  v 


in  o  n  ro  cn 


O 


CD 

0) 

-O 


c\j  iv.  c\j  r-  cm 
r***  co  CO  «—■  rn 


o  w  n  in 
-  rvj  ru  m  cn 


(D 

> 


1  I  t  I  1 


(N.  W  N  (Vi  N 
m  N  CD  to  -* 


’TCDC'JN  — 

—  — »  cvj  cm  m 
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APPENDIX  C 
SOFTWARE  DIRECTORY 

This  Appendix  gives  a  listing  of  the  various  programs  used  in  this  thesis.  Each  program 
listing  contains  a  summary  of  how  the  program  is  used,  user  inputs,  program  outputs  and 
any  additional  features.  Programs  are  listed  in  two  groups;  Nusselt  number  distributions 
and  exit  mixed-mean  temperature  measurement.  All  programs  were  written  in  BASIC  4.0 
for  use  on  the  HP9836S  computer. 

I.  NUSSELT  NUMBER  DISTRIBUTION  PROGRAMS 

A.  TCHECK: 

This  program  acquires  multiple  channel  thermocouple  data  and  performs 
temperature  calibration  equations  to  provide  measurements  in  Degrees  Celcius. 
user  inputs:  thermocouple  voltage  readings,  X 

program  output:  provides  monitor  summary  of  thermocouple  number,  data  aquisition 
channel  number,  voltage  reading  and  corresponding  temperature 
reading  in  Deg.  C. 

B .  TEMCURV: 

This  program  aquires  multiple  channel  thermoucouple  data  and  creates  two  files 
to  be  used  in  other  programs  for  the  curved  heating  channel.  Primary  purpose  is  to  supply 
data  for  Nusselt  number  calculations, 
user  inputs:  a)  thermocouple  voltage  readings. 

b)  run  number. 

c)  voltages  across  shunts  for  all  four  heaters,  mvolts. 

d)  voltages  across  heaters,  volts. 
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e)  ambient  pressure,  inches  of  mercury. 

f)  pressure  drop  across  orifice  plate,  inches  of  water, 
program  outputs:  a)  data  file,  Idata,  provides  data  required  to  determine  heat 

input,  conduction  loss  and  Dean  number;  Runno,  Ampl,  Amp2, 
Amp3,  Amp4,  Vhl,  Vh2,  Vh3,  Vh4,  Pamb,  Deltap,  Tmi, 
Tccamb,  Tcvamb,  Ro. 

b)  data  file,  Tdata,  provides  temperature  data  in  Deg.  C  for 
determination  of  local  Nusselt  numbers. 

C.  NUCURV3B: 

This  program  determines  heat  transfer  coefficients  and  Nusselt  numbers  by 
conducting  surface  energy  balance.  Heat  flux  inputs  and  conducted  heat  losses  are 
determined. 

user  inputs:  a)  reads  data  files  Tdata  and  Idata  from  program  TEMCURV. 

b)  estimated  Dean  number. 

program  outputs:  a)  data  file,  Hdata,  of  heat  transfer  coefficients,  local  Nusselt 

numbers  streamwise  and  span  wise  locations,  measured  and 
corrected  temperatures  as  follows;  H,  Nus,  X2,  Z,  T1  and  T. 
additional  features:  a)  provides  spanwise-averaged  Nusselt  numbers. 

b)  using  energy  balance,  calculates  local  mean  temperatures. 

c)  with  estimated  Dean  number,  performs  iteration  process  to 
determine  actual  Dean  number. 

D.  NUSCCPL2 

This  program  reads  data  file  and  plots  spanwise  variations  of  Nusselt  numbers 
by  rows  for  concave  section  of  curved  heating  channel  utilizing  thermocouples  31-115. 
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user  inputs:  local  Nusselt  numbers  along  with  corresponding  stream  wise  and  spanwise 
distances  (output  data  file,  Hdata,  from  program  NUCURV3B). 
program  outputs:  Nusselt  number  distribution  contours  for  all  rows  of  concave  curved 
section. 

E.  NUSCVPL2 

This  program  reads  data  file  and  plots  spanwise  variations  of  Nusselt  numbers 
by  rows  for  convex  section  of  curved  heating  channel  utilizing  thermocouples  116-120. 
user  inputs:  local  Nusselt  numbers  along  with  corresponding  streamwise  and  spanwise 
distances  (output  data  file,  Hdata,  from  program  NUCURV3B). 
program  outputs:  Nusselt  number  distribution  contours  for  all  rows  of  convex  curved 
section. 

II.  EXIT  MIXED-MEAN  TEMPERATURE  MEASUREMENT 
A .  TEMTRAV 

Program  for  control  of  temperature  traverse  probe  positioning  and  data 
acquisition  in  the  curved  heating  channel.  Interfaces  with  MITAS  controller  to  control  the 
probe  position.  Provides  temperature  difference  between  exit  temperature  and  ambient  inlet 
temperature. 

user  inputs:  a)  estimate  of  Dean  number 

b)  initial  position  of  pressure  probe. 

c)  number  of  y  and  z  points  to  take  data  at. 

d)  resolution  of  data. 

e)  ambient  inlet  temperature  channel  number. 

f)  temperature  probe  channel  number. 

g)  ambient  pressure. 
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h)  pressure  drop  across  orifice  plate. 

i)  delay  time  to  start  the  experiment 

program  outputs:  a)  data  file.  Temp,  provides  position  and  temperature  data  as 

follows:  Y,  Z,  temperature,  inlet  temperature,  temperature 
difference. 

b)  provides  actual  Dean  number  through  iteration  process. 

B .  TMXDMPLOTB 

Program  to  plot  time  averaged  traverses  of  temperature  probe  measurements. 
Will  plot  up  to  ten  radial  traverses  for  different  spanwise  locations, 
user  input:  a)  temperature  data  file,  TEMP  (output  from  program  TEMTRAV). 

b)  spanwise  positions  to  plot 
program  output:  temperature  traverse  plot 

C.  TMXDMEAN 

Modified  version  of  TMXDMEAN  to  plot  traverse  of  mean  temperature,  mean 
velocity  products.  Will  plot  up  to  ten  radial  traverses,  for  different  spanwise  locations. 
Program  has  been  modified  to  determine  exit  mixed-mean  temperature  using  process 
described  in  Chapter  3,  Section  C,  part  2. 

user  inputs:  a)  temperature  data  file,  TEMP  (output  from  program  TEMTRAV). 

b)  velocity  data  file,  VEL  (output  from  program  VELOCITY). 

c)  spanwise  positions  to  plot. 

program  output:  a)  mean  temperature,  mean  velocity  traverse  plot 
b)  monitor  output  of  exit  mixed-mean  temperature. 
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DATA  FILE  LISTING 


A  Summary  of  the  data  files  and  corresponding  experimental  parameters  is  given  in 

two  parts,  the  first  of  which  involves  the  determination  of  all  Nusselt  number  results 
and  the  second  with  the  calculation  of  the  measured  exit  mixed-mean  temperature 
results. 

A.  NUSSELT  NUMBER  RESULTS. 

The  following  lists  of  data  were  used  for  the  determination  of  all  Nusselt  number  data 
presented  in  this  thesis.  This  includes  both  the  spanwise-ave raged  Nusselt  data  and  the 
local  spanwise  Nusselt  number  distributions. 


experimental  parameters: 

Dean  nnumber  58.9 

74.2 

98.6 

126.4 

150.9 

174.4 

198.9 

4 

Pressure  drop 

.031 

0.05 

0.09 

0.15 

0.215 

0.29 

0.38 

(in  water) 

Heat  flux/heater 

38.8 

37.8 

53.75 

52.9 

73.6 

72.2 

76.0 

(watts) 

Max.  surface  temp. 

56.5 

48.6 

54.9 

49.3 

57.1 

51.5 

50.2 

(  TC#199,  OC  ) 

data  file  listings; 

Tdata/ldata  # 

22-24 

26-28 

6-8 

9-11 

1  -3 

4-6 

7-9 

File  number 

9 

9 

9 

9 

10 

10 

10 

Hdata  # 

50 

75 

100 

125 

150 

175 

200 

(file  #19) 


* 
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B.  MEASURED  MIXED-MEAN  TEMPERATURE  RESULTS 


The  following  list  of  data  provides  the  data  information  used  in  calculating  the 
measured  exit  mixed-mean  temperature  and  the  plotting  of  the  temperature  traverse 
contours.  Also  included  is  the  velocity  file  number,  which  is  presented  by  Hughes  [Ref. 
8]  and  Tuzzolo  [Ref.  14]. 

experimental  parameters: 


Dean  number 

58.5 

73.7 

98.1 

126.0 

149.8 

173.5 

198.4 

Pressure  drop 
(  in  water) 

0.031 

0.05 

0.09 

0.150 

0.215 

0.29 

0.38 

Heat  flux/heater 
(watts) 

40.0 

47.2 

56.7 

63.3 

73.2 

81.9 

92.2 

Max.  surface  Temp. 
TC  #199,  oc  ) 

60.4 

60.3 

59.2 

56.8 

57.6 

58.8 

59.7 

data  file  listing: 

Tdata/ldata  # 

3.4 

5.6 

7,8 

9,1  0 

11,12 

0,1 

2,3 

File  number 

1  5 

1  5 

1  5 

1  5 

1  5 

1  7 

1  7 

Temp  # 

50 

75 

100 

125 

150 

1  75 

200 

File  number 

1  5 

1  5 

15 

1  5 

1  5 

17 

1  7 

Vel  # 

1 

1 

3 

4 

4 

5 

51 
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